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ABSTRACT 

This is the first paper about the fragmentation and mass outflow in the molecular cloud by 
using three-dimensional MHD nested-grid simulations. The binary star formation process is 
studied paying particular attention to the fragmentation of a rotating magnetized molecular 
cloud. We assume an isothermal rotating and magnetized cylindrical cloud in hydrostatic bal- 
ance. Non-axisymmetric as well as axisymmetric perturbations are added to the initial state 
and the subsequent evolutions are studied. The evolution is characterized by three parameters: 
the amplitude of the non-axisymmetric perturbations, the rotation speed, and the magnetic 
field strength. As a result, it is found that non-axisymmetry hardly evolves in the early phase, 
but begins to grow after the gas contracts and forms a thin disk. Disk formation is strongly 
promoted by the rotation speed and the magnetic field strength. There are two types of frag- 
mentation: fragmentation from a ring and that from a bar Thin adiabatic cores fragment if a 
thickness is smaller than 1/4 of the radius. For the fragments to survive, they should be formed 
in a heavily elongated barred core or a flat round disk. In the models showing fragmentation, 
outflows from respective fragments are found as well as those driven by the rotating bar or the 
disk. 



Key words: binaries: general - 
stars: formation. 



1 INTRODUCTION 

It is important to understand how a binary star system forms from 
a molecular gas cloud because most main-se quence stars are in 
binary or multiple systems (e.g. lAbt„,1983.: .Duguennov & Mavod 
I199II) . Furthermore, recent observations indicate that compan- 
ions may be even more common for pre-main sequence stars 
iRichichi et al. 1994). Numerical simulations show that it is diffi- 
cult for single stars to capture an other one and to form a b inary sys- 
tem in the main sequence phase jKroupa & Burkeri200lh . Accord- 
ingly, most stars seem to be born as binaries rather than single stars 
in a gas cloud, and molecular bipolar outflows are often observed 
in star forming regions. It has been shown using two-dimensional 
magnetohydrodynamical simulations that outflow plays an impor- 
tant role in the star formation process because the excess angular 
momentum, which is conserved fro m the host cloud, is rem oved by 
these outflows from the formed core iTomisaka 2000, 2002). These 
outflows are caused by the twisted magnetic field lines made by the 
rotation near the adiabatic core. In order to investigate the cloud 
fragmentation and later binary star formation, three-dimensional 
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calculation is required. In this letter, we discuss fragmentation in 
the course of star formation which leads to binary star formation. 

Fragmentation of isothermal cloud is studied by many au- 
thors usin g three-d imensional hydrodynamical simulations (e.g. 
feodenhei mer et afl ^00). The criterion for an initially uni- 
form cloud to fragment is described by ao (thermal-to- 
gravitational energy ratio) of less than ~ 0.2 — 0.5 if we re- 
strict ourselves to the isothermal gas. The clou d with qq > 
0.2 — 0.5 experiences the run-away collapse iLarsonI 1 19691 
iNorma n. Wilson & Barton 1980). When the gas density increases 
up to n > 10^" cm ~^, the isothermal gas approx i matio n 
breaks down. Recently iBate. Bonnell & Bromri] |200l l2003h . 
ICha & WhitworthH2003l) and lMatsumoto & Hanaw j <2003bh have 
studied the cloud evolution from isothermal to adiabatic state us- 
ing a barotropic equation of state and they clarified the evolution 
of the cloud with no magnetic field. However, there are few papers 
on the fragmentation of the magnetized cloud. Ibos3 i2002l) studied 
the fragmentation taking account of radiative transfer with the Ed- 
dington approximation instead of barotropic equation of state, but 
simplified the effect of magnetic field. In his studies, the angular 
momentum transfer by the magnetic field is ignored. However, the 
angular momentum is largely removed by magnetic braking and 
molecular outflow in the cloud evolution process iTomisaka 2000). 
Therefore, to explore the fragmentation of the magnetized cloud. 
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full three-dimensional magnetohydrodynamical MHD simulations 
are necessary. 

Here, we employ the nested grid MHD code, which always 
maintains sufficient spatial resolution in the central region. The 
nested grid method is powerful to study the star formation process, 
because this requires a large dynamic range in spatial dimensions. 
In the simulations, we use the barotropic equation of state instead 
of inclusion of radiative transfer. In this study, we calculate the non- 
axisymmetric dynamical contraction of the cloud from 4 x 10^ 
cm^'^ to 1 X 10^^ cm^'' in number density, and investigate how 
fragmentation proceeds and outflow occurs in the course of binary 
star formation. 



2 MODEL AND NUMERICAL METHOD 

To study the star formati on (binary star forma tion and their outflow) 
in the filamentary cloud ( Mizu no et al. , we consider a cylin- 

drical isothermal cloud in hydrostatic balance as the initial condi- 
tion. The cloud is assumed to rotate parallel with a rotation axis 
which coincides with the cylindrical axis and the magnetic field 
lines are assumed to run parallel to the rotation ax is, because mag- 
netic field line s often turn to the cylindrical axis iGoodman et al. I 
[l990.: Ward-Thompson et al. 2000). Using the cylindrical coordi- 
nates (r, (j), z), the density, azimuthal velocity, and magnetic flux 
density distributions in the radial direction are taken to be as fol- 
lows ( Stodolkiewicz 1963): 

po{r) = p,[l + {r^/8H^)]~\ (1) 

v^o{r) = rn,[l + {r'/8H^)]''^^\ (2) 

B.oW = B,[l + {ry8H^)]~\ (3) 

where pc, fie, and Be represent the density, angular rotation 
speed, and magnetic flux density at the center of the cylindrical 
cloud, respectively, and H is the scale-height as — (c^ + 
B'^/8iTpc)/{^nGpc — 2f2^)i where Cs denotes the isothermal 
sound speed. Molecular gas obeys the isothermal equation of state 
below Peri ~ 10^° cm~'^ but it becomes adiabatic if p > Pcri 
iToh line 1982). To mimic this, we adopt a two-components equa- 
tion ofstate_as_P_^^_c^^_J-c^ iBonnell & Bata 
Il994 iBate. Bonnell & Brommll20o3) . The initial central density 
and the critical density are chosen as pco = 4 x 10^ cm~^ and 
Pcri = 10'"cm-^ 

To initiate contraction, we added axisymmetric and non- 
axisymmetric perturbations to the above magnetohydrostatic equi- 



librium as 

p(r,z,0) = poir)[l + Sp4z)][l + Sp4,{r,4>)], (4) 

B,(r,z,(P) = B,o{r,<l,)[l + SB4,{r,<f>)], (5) 

with 

Spziz) = A;, cos{2tv z/L), (6) 

^ , ^_ , ( Asir / H)^ cos(md>), iorr^H^^ 



where L is chosen equal to the wavelength of the most 
unstable Jeans mode acquired by the linear analysis 
^ Matsumoto, Nakamura & Hanawa 1994). 

In this letter, we restrict ourselves to the relative amplitude of 
the axisymmetric perturbation = 0.1. As the non-axisymmetric 
perturbation, only m=2 mode is included. Equations J4} and js} 
ensure the ratio of the density to the magnetic flux constant in 



the azimuthal direction. Models are parameterized with three non- 
dimensional parameters: the amplitude of the non-axisymmetric 
perturbation as A^p, the magnetic-to-thermal pressure ratio as a(= 
B^/A-Kpc(?s) and the angular speed normalized by the free-fall 
timescale as lj = n^/i'i-KGpaf''^. We calculated 51 different 
models with A,^ = (0, 0.01, 0.1, 0.2), q = (0, 0.01, 0.1, 1, 5), 
andoj = (0, 0.1, 0.5, 0.7). 

In order to solve the central region with higher spatial resolu- 
tions, we adopt an MHD neste d grid method b ased on the Carte- 
sian coordinate (for detail, see iMachiddEoOSI) . In this method, a 
number of grids with different spacings are prepared, in which 
finer grids cover the central high-density portion and the coarser 
ones cover the cloud as a whole. Each grid has successively dif- 
ferent cell widths by a factor two. We use 128 x 128 x 32 cubic 
cells in the x-, y- and z-directions. We added a new finer grid to 
maintain the Jeans condition, Aj/4 > h, with an ample margin 
(Truelove et al. 1997), where Aj and h are the Jeans length and 
the cell width, respectively. Whenever one-eighth of the minimum 
Jeans length becomes smaller than the cell width of the finest grid, 
a new finer grid is added to the nested grid. Total number of cells 
is 128 x 128 x 32 X ^max, where ^niax(^ 17) denotes the maxi- 
mum grid level. Only when a high density fragment escapes from 
the region covered by the finest grid, the Jeans condition is violated 
in our simulations. 

The MHD code u ses the numerical fluxes proposed by 
iFukuda & Hanawj J1999') but modified to solve the isothermal and 
polytropic gas. We adopt the MUSCL approach and predictor- 
corrector method in the time integration. The Poisson equation 
is solved by the multigrid iteration method on the nested grid 
iMatstmioto & H anawa 2003a). We adopt mirror symmetry with 
respect to the z = and z = L/2 planes. We also set the fixed 
boundary condition on the surface of r = 2H. 



3 RESULTS 

The cylindrical gas cloud collapses to become spherical in shape 
in the early collapse phase. If the initial gas cloud has neither 
magnetic field nor rotation speed, the cloud continues to collapse 
spherically a nd it forms in to a small massive core in the center 
of the cloud jLarsorJI 19691) even in a cylindrical cloud jTomisak^ 
Il995t) . In the case of a gas cloud with magnetic fields or rota- 
tion, it collapses to form a pseudo-disk due to the magnetic pres- 
sure or centrifugal force ( Tomisaka' ' 1 995*; 'Nakamura & Hanawal 
.1997 ; Matsumoto. Hanawa & Nakamura 1 997). If only the axisym- 
metric perturbation exists, a round disk is formed, but it does 
not fragment as long as the axisymmetry holds. When the non- 
axisymmetric perturbation is added, the cloud evolves to form a 
non-axisymmetric shape, and fragments later for some range of 
the parameters. In any case, the central density (pc) increases with 
time and the gas becomes adiabatic when pc > W^'^ cm""*. For 
convenience, we divide the evolution into two phases: collapse 
phase (pc < lO^'^cm^'^) and accretion phase {pc > 10^'^ cm^'^). 
In this equation of state, the temperature gradually increases near 
p < Pcri, however, the solution acquired in p < pcri is almost the 
same as that obtained for the isother mal collapse. The outfl ow is 
driven by the twisted magnetic fields iTomisak J 1 998l l2002h after 
~ 10'' yr have passed from the core formation epoch pc = Pcri- 

Fig. Q shows the pole-on views of the adiabatic cores (upper 
panels) and the side-views of the outflow regions (lower panels) 
for three models in the accretion phase. The thick line denotes the 
boundary of the adiabatic core (p ^ pcri) for the upper panels. 
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In the lower panels, the thick lines which represent the isoveloc- 
ity lines of — 0.1, 1 and 4 kms^^ indicate the outflow re- 
gion. In Fig.^, a model without non-axisymmetric perturbations 
(A^j, = 0) is shown. Model parameters {a,u)) — (0.1,0.1) cor- 
respond to the cloud with a weak magnetic field and slow rotation. 
This cloud collapses almost spherically and forms a small spherical 
core at the center of the cloud. We have no evidence of frag menta- 
tio n in thi s model. On the other hand, in the studies of Boss ( 2002') 
and lLi & 'Nak amura ( 2002), fragmentation occurs in this parameter. 
However, their results are not inconsistent to ours, taking account 
of the fact that the magnetic braking is ignored in their calculation. 
We also found a strong outflow is driven near the adiabatic core 
after the central density reaches ~ 10^'^ cm~"^ (lower panel). 

A model with small non-axisymmetric perturbation (A^ = 
0.01) is shown in the middle panel (h):(a,ui) — (0.01,0.5). In 
this model, a pseudo-disk is formed in the isothermal collapse 
phase, because the gas cloud is rotating rapidly (u = 0.5). The 
disk evolves into a ring in the accretion phase, and fragmenta- 
tion occurs on the inner edge of the ring. Two contracting cores at 
{x, y) = (75 AU, 125 AU) and (-75 AU, -125 AU) are formed. 
The right panel (c):{a,u)) — (1,0.5) shows the model with a 
large amplitude of non-axisymmetric perturbations (A^ — 0.2). 
In this case, non-axisymmetric patterns grow sufficiently in the 
isothermal collapse phase and a bar is formed at the beginning of 
the accretion phase. In the accretion phase, the bar grows longer 
and thinner. The bar fragments and forms two contracting cores 
lix,y)={~10 AU, -25 AU), (10 AU, 25 AU)]. Fig.Qshows that 
there are at least two kinds of fragmentation: fragmentation from a 
ring (b) and from a bar (c). 

Fig. Q;- shows that lobes of outflow are ejected from each 
contracting core {y ~ ±25AU), whose velocities reach ~8 
km s~^. As for the outflow, low-velocity outflows are seen also 
in Fig. each of which is connected to the contracting core. 
In both the bar and ring fragmentations, two types of outflows 
co-exist: the inner one is connected to the cores formed by the 
fragmentation and driven by their spin motion, and the other is 
driven by the rotation motion of the ring or the bar. However, the 
strength of the outflow is different for respective fragmentation 
patterns. Fast outflows appear in the bar fragmentation (c), while 
weak ones appear in the ring fragmentation (b). This difference 
comes from the angular momentum distribution between the spin 
motion of the cores and their orbital angular momentum. In the 
ring fragmentation, the fragments have the angular momentum in 
the form of orbital angular momentum. On the other hand, in the 
bar fragmentation, the angular momentum is almost evenly shared 
between the spin and the orbital angular momenta. In the model 
of Fig.Q;, the fragments are expected to merge each other due to 
their small orbital angular momentum. 



4 DISCUSSION 

As mentioned in §|3| modes of fragmentation and shapes of the 
outflow depend on the growth of non-axisymmetry. This shows that 
the cloud evolution can be classified by this growth. In order to 
characterize the shape of the central region, we define oblateness 
(Sob) and axis ratio (Sar) as follows: 

eob = \/ hi ■ hs/hz, (8) 
Ear = hi/hs, (9) 



where hi, ha and h^ mean the length of the long-, short-, 
and g-axis for the gas wit h p > (l/10)pc, respectively (see 
iMatsum oto & Hanawa''l999. for definition of h). 

In Fig.|2| the oblateness and the axis ratio of the central re- 
gion are plotted for all models at the beginning of the accretion 
phase (pc = Peri = 10^" cm~'^). Curves in this figure illustrate 
the evolutionary tracks for four different models (A^, a, u!)= (0.01, 
0.01, 0.5), (0.2, 0.1, 0.5), (0.2, 1, 0.5), and (0.2, 5, 0.5). The clouds 
evolve from the lower-left corner, which means the early shape is 
almost spherical, to the upper-right region via the upper-left region. 
First, the clouds evolve vertically upward. This indicates that the 
oblateness increases. Then, the curves change their directions to 
the right for three models with A^ = 0.2. Evolutionary tracks of 
the high-density region in Fig. |2| show very clearly that the non- 
axisymmetric perturbation grows only after the disk is formed. That 
is, the axis ratio increases after the oblateness grows considerably 
large (Sob ^ 4). Because the gas is partially supported in the radial 
direction by the magnetic pressure and the centrifugal force, the ra- 
dial contraction is more delayed than that of the ^-direction. As a 
result, the disk is formed (Sob > 4) earlier for models with larger 
Q and uj. In such models, the non-axisymmetry grows sufficiently 
and forms a large bar in the isothermal collapse phase. However, 
the non-axisymmetry hardly evolves in the model with small a and 
Lu, because it takes a long time for disk formation. 

Adiabatic cores are classified into three types: core (Fig^: 
Eob < 4 and £ar < 4), disk (FigQJ): £ob > 4 and £ar < 4), 
and bar (FigQ;: Eob > 4 and Ear > 4). Some disks evolve into 
the rings, in which the central density is lower than that outside, 
even if it has the same oblateness and the axis ratio. Models are 
divided into three kinds of evolution: thin-line symbols mean that 
the models form cores which continue to collapse, while thick-line 
and filled symbols represent the models in which fragmentation oc- 
curs. Although the fragments survive against mutual merger in the 
models of thick-line symbols, they merge each other in the models 
of filled symbols. 

Since the models with thick-line or filled symbol are dis- 
tributed only in the region Eob > 4, it is concluded that the nec- 
essary condition for fragmentation is Eob > 4, which means that 
a thin disk or bar is necessary for fragmentation. Even in the re- 
gion of Eob > 4, there are exceptional three models in which the 
adiabatic core evolves to form an almost axisymmetric thin disk 
or ring because non-axisymmetry is counteracted in the isothermal 
collapse pha se. The shapes of the bar a nd ring are also seen in 
IBossI <2002') and lLi&Nakamurj ilOOft. It is concluded that for 
the adiabatic core to fragment, the central gas needs to form a suffi- 
ciently thin disk or bar in the isothermal collapse phase. The thick- 
line symbols are distributed in two regions enclosed by circles. The 
fragmented cores in the circles seem to evolve into a binary star be- 
cause the length of the semi-major axis expands or oscillates. Mod- 
els in which the fragments merge with each other are distributed 
outside the domain. Outside of the domain, a few fragments sur- 
vive at the end of the calculation but the separation between the 
two cores shrinks with time. These will not grow into binary stars. 

In conclusion, if the following criteria are fulfilled at the end of 
the isothermal collapse phase, fragmentation occurs and fragments 
survive to form binary or multiple stars. 

(i) Eob > 4 : fragmentation condition. 

(ii) < Ear < 2, or Ear > 10 : survival condition against 
merger. 

The symbol in this figure indicates the shape of the adiabatic core 
when the fragmentation occurs or the calculation ends. Fig.|2|con- 
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firms the idea that the fragmentation occurs from a ring (Q) or bar 
(A) and not from a core (□). 

We showed the domain which may evolve into binary stars 
in Fig. |2| Outside the binary-forming domain, the core formation 
model (□) results in a compact massive core at the center of cloud. 
The bar-fragmentation model (A) in which the mutual merger of 
the fragments occurred also leads to a disk-like core at the center 
of cloud. This configuration seems to lead to single star formation. 
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Figure 1. Three typical models of the adiabatic core are plotted (core- 
type: left, disk- or ring-type: middle, bar-type: right). Model parameters are 
{A^, a, uj) = (a): (0.0, 0.1, 0.1) , (b):(0.01, 0.01, 0.5), and (c):(0.2, 1.0, 0.5). 
Density (contour and false color) and velocity vectors (arrows) are plotted 
on a plane including the major axis of the core and the z-axis (lower pan- 
els) and on the z=0 plane (upper panels). The elapsed time, the maximum 
density, grid level and scale length are also plotted. Thick lines in the up- 
per panels denote the contour line of p = lO'^'^cm"^ showing the outline 
of the adiabatic core, while those in the lower panels mean the isovelocity 
curves representing the outflow region (v^ = 0.1, 1, and 4 km s"^) 
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Figure 2. The axis ratio and oblateness are plotted when the gas at the 
cloud center becomes adiabatic pc = Pcri- The abscissa indicates the axis 
ratio and the ordinate is the oblateness. Both are calculated for the gas 
with p > (l/10)pc. Fragmentation occurs above the horizontal dashed 
line (Eob > 4). The domains in large circles mean that fragmented cores 
survive to form binary or multiple stars. The solid lines represent the evo- 
lutional tracks for some typical parameters. Dotted lines represent the evo- 
lution after the adiabatic core is formed. Asterisks show the points of frag- 
mentation. The symbols □ (core), "O" (disk), Q (ring) and A (bar) denote 
the shape of the adiabatic core when the fragmentation occurs or the cal- 
culation ends. Thick-line and filled symbols represent the cores which ex- 
perience fragmentation. Thin-line symbols represent the non-fragmentation 
models. However, filled symbols indicate merger core after fragmentation. 



